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A quark confinement potential is derived from QCD Lagrangian under the two assumptions in the 
energy region relevant to quark confinement. The possibility of perturbation in QCD is assumed, 
and the momentum of the gluon transfered by the quark is larger than the mass of the internal 
quark, so that the expansion of m 2 /q 2 is assumed to be reliable. This potential comes from the 
one loop diagrams in the expansion of m 2 /q 2 , especially here, the vacuum polarization diagram is 
examined. This explains the quark confinement mechanism on the basis of QCD and the Regge 
trajectories which is governed by a linear potential. 
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Quarks are known as the constituents of the nuclcon and the meson, and their interaction is believed to be explained 
by quantum chromodynamics (QCD). However there is no direct experimental evidence of detecting an isolated single 
quark. So it is said that quarks arc always confined in the hadron. QCD tells us that quarks behave as weakly 
interacting particles in high energy states, which is referred to as "asymptotic freedom". In this energy states, 
perturbative calculations are possible so that perturbative QCD has been actively researched Q. In order to motivate 
\ this letter, it is useful to cite "If indeed the quarks behave approximately like free particles and their masses are rather 

■ small, the critical question, then, is why don't we see free quarks in the final state? This is the well-known problem 
I/") [ of quark confinement." In low energy physics, however, perturbation is questionable because of the behaver of 

■ the larger coupling constant of QCD. So quark degree of freedom doesn't appear manifest. It is possible to describe 
nuclear phenomena with the nucleon and the meson, i.e., meson theories instead of quarks. In order to describe low 
energy nuclear phenomena with quarks, several potential models are used such as bag models, harmonic oscillator 
potential, linear confinement potential model and so on |}|^| . Another kind of model is the string model from which 
a linear confinement potential is derived naturally, but in which the concept of a particle should be changed into a 
string Q. All those potential models have their own quark confinement conditions. 

If QCD is a true theory for the interaction of quarks and gluons, why isn't the quark confinement phenomenon 
— — I | explained well with it? The hint to this question is the Lamb shift. With this hint, it is shown that QCD has a certain 
confinement potential which can be derived from QCD Lagrangian only. Consequently, this potential gives the reason 
why quarks are always confined in the hadron. Quark confinement is explained by the one loop effect of QCD, as if 
the Lamb shift is done by that of quantum electrodynamics (QED). The essential difference between QCD and QED 
in the two phenomena is the magnitude of the coupling constants rather than the difference between the Abelian and 
non-Abelian gauge theories. In QED, a hydrogen atom is the well-known bound state of a proton and an electron. 
The coupling constant of QED is so small compared to that of QCD that the size of the hydrogen atom should be 
much larger than that of the nucleon which is a bound state of quarks. The larger coupling constant will cause the 
larger momentum of the gluon which is transfered by the quark in order to localize quarks in the size of a nuclcon. 

In the calculation of the Lamb shift, the momentum of the photon transfered by the electron is so small that the 
one loop diagrams are calculated in the expansion of q 2 /m 2 , and the result is the famous splitting of the 2p and 2s 
states of the energy levels of a hydrogen atom. This splitting is generated by the following effective potential (|||: 

a 4a 2 

V(r) =--- -j^S(r) (1) 

in configuration space. In momentum space, this potential corresponds to the propagator of the photon which is 
modified by the one loop effect, i.e., the vacuum polarization: 

_ (g^ - q/y) e 2 q 2 

^~ q 2 { m^m 2 '- [2) 



*e.mail : yongkyu@phya.yonsei.ac.kr 



1 



However what is calculated in such an expansion in QCD may be inappropriate to explain the hadronic phenomena, 
because quarks are interacting so strongly with each other. The quarks confined in the hadron may not be in low 
energy-momentum states because of the large coupling constant which will produce a large momentum transfer. So the 
basic assumptions of this letter are that the quarks inside a hadron are so relativistic that perturbation is meaningful 
and the expansion of the inverse momentum of the gluon, that is, m 2 j q 2 is reasonable rather than q 2 /m 2 . 
The QCD Lagrangian is 



FZ» = MS - d.A^ + gf abc AlAl , (3) 

where the index a runs 1,2,.. ..,8. The effective potential to one loop order can be calculated from the Feynman 
diagrams of Fig. [IJ Those diagrams contribute to the following transition amplitude: 

M = ff 2 u(p 2 ){ 7 ^ 1 (g 2 ) + £-F 2 (q 2 )a^ q a }T a u( Pl )D^Jt, (4) 

where J h v means the quark current which consists of the bispinor of an antiparticle for the meson or the composite 
of the bispinors of two quarks and gluons for the nucleon. This amplitude is not a direct observable in the case of 
a bound state but an analogy of the amplitude of a scattering process. As the purpose of this letter is to show a 
qualitative confinement potential form QCD, the vacuum polarization diagram is considered only. It is calculated as 



n«l(q 2 ) = ig 2 ^- d Tr(T a T b ) J ^~ 



d d k 1 1 



q 2 nr5„h , o . r l 7 1, ,— m 2 , 

= —vk^ {q g ^ ~ q ^ ){ 7 2 2 l0g W )} + finitc(9 } ' ( } 

where rif is the number of flavors to contribute to the loop of the quark The infinite part in the above result 
contributes to the renormalization of the wave function of the gluon. The remaining finite part is calculated as 

TT ( 2s 9 2n fSab (2 , f 1 , n m f x{l-x)q 2 -m 2 
Ilfinito {q ) = -~^y— (q 9^v - q^qv) / dxx(l - X)\0g{ g ) 

4-7r z ,/n — m z 



If it is expanded in —Am 2 /q 2 which is called the expansion of m 2 /q 2 throughout this letter, it is easy to calculate the 
effective potentials. This vacuum polarization modifies the gluon propagator as 

nab( 2s _ 5 ab {9^ - g'V) S ac (g^ - g»g») d 5 bd (g™ - q° g") 

lLv\q ) 2 2 per 2 

= q 2 [1 + 1^ 18 - 6 l0g( -^ } - V 

)2 + ( f )2l ° g( -f ) + } + 1' (7) 

where the coupling constant g is renormalized one from now on, since the infinite part of the vacuum polarization is 
canceled by the suitable counter term which comes from the redefinition of the physical parameters, i.e., wavefunctions, 
m,g and so on. In the center of mass frame of a quark and the remaining part, the four momentum of the gluon 
transfered by the quark is reduced to a three momentum, because the energy of the quark before the interaction with 
the gluon is exactly equal to that after the interaction. Now the Fourier transforms of each term are 



(27r) 3 q 2 47rr' 

/ 7^4log(qV q - r --^log( 7 r), 
J {Ztt)- 1 ' q z Zirr 
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1 V\ 

log(q 2 )e*»" - — {- - log( 7 r)}r 3 , (8) 

where the table of integrals of Ref. |tJ is used and 7 is the Euler constant in the second Fourier transform. The first 
one is the usual Coulomb potential as encountered in QED. The second one is known as the Uehling potential [||. 
The third one is the linear potential which confines quarks in the hadron forever. The last two equations are only for 
reference, which give higher order confinement potentials. Alternatively, the linear confinement potential can be also 
estimated from the exact equation (Q) by ignoring the last term which behaves approximately like log(q 2 ) at large 
momentum in the gluon propagator, so one can check that the integration interval in the Fourier transformations and 
higher power terms in Eq. (Q) do not matter. From the Fourier transform of the propagator of the gluon, the effective 
potential by which a quark is influenced is calculated as 

V(r) = - Ci A + + ^^log( 7 mr) + °-fl m W s r + , (9) 

r lair r 07T r 2ir 

where cf is the color factor which is | for the meson and | for the baryon and a s — j- The non- Abelian nature of 
QCD gives only a constant, such as the color factor from the trace in Eq. The diagrams Fig. f-(e) and (f) which 
are quite different from QED do not give any confinement potential in the gluon propagator but give only Uehling 
potentials because there is no internal quark line. The complete effective potential may be needed to consider further 
through out the remaining diagrams of Fig. f and bremsstrahlung diagrams. 

There is the string model that has such a linear potential Q . This model explains a linear dependence of the angular 
momentum of the states of highest angular momentum, on the square of the mass, namely, J = a'E 2 + constant. 
This relation holds for the case of a constant energy density k of a rotating string, that is, to a potential of the form 
V = kr. The value of k is obtained from the slope a' with the relation a' = l/2nk, which is derived in Ref. [0. 
Inserting the observed value a' — 0.93 GeV 2 from the plot of spin J against mass squared for baryon resonances of 
the A family ( s = 0, I = | ) and the A family ( s = -1, I = ), the coefficient k is 0.85 GeV fm _1 . Comparing it 

with Eq. (||), then the coefficient k is matched to the result of this work as k = c ^" a S/=u d a m /' w here the mass 
squared is summed up to the quark satisfying the condition —Arnj/q 2 < 1. Therefore the Regge trajectories can be 
regarded as a strong experimental evidence of the linear confinement potential. 

There is a theory that has a large coupling constant but not a confinement potential. It is the very N-N potential in 
nuclear physics which is produced by the force mediated by pions, vector mesons and so on by coupling to nucleons. 
The coupling constant of g^ NN /(A%) is about 14.3 in the literature jpTof] . The pion mass is the reason why there is no 
such a potential. The pion propagator corrected by the one loop effect in the expansion of m 2 /q 2 is 



q 2 — m 2 q 2 — m 2 q 2 — m 2 



1 m , -9 : 
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- 2 {l + f^{C 1 + C 2 log(q 2 ) + C 3 — + } + ], (10) 

q A — m z 47r z q z 

where the free propagator should not be expanded because of renormalizability. The Fourier transform of the term 
which has the coefficient C3 and is expected to give a confinement potential like Eq. (0) still has an exponential decay 
factor due to the mass of the pion as followings: 



d 3 q m 2 4q . r _ J_ 



(27r) 3 q 2 (q 2 + m 2 ) 4nr 



(l-e-" r ). (If) 



The higher order terms of m 2 j q 2 also do not give a confinement potential, which can be proved by differentiating the 
above Fourier transform twice with respect to r. Therefore the N-N potential has no confinement potential, though 
it has a large coupling constant. 

However there remains still a question: how large coupling constant produces such a confining potential? Since such 
a potential can be also derived in QED, if the possibility of the expansion of m 2 /q 2 is forced to assume, it is necessary 
to justify the assumptions of this letter. The running coupling constant of QCD is derived in the perturbative region, 
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so in the non-perturbative region, it is useful to use that the velocity of the electron v is equal to the coupling constant 
a in the Bohr model. Using the relativistic energy momentum relation p 2 H = Mfj — (p e + p p ) 2 in the center of mass 
frame, the identity is calculated as 



In I (m e + m v ) x \(M'h — m 2 — m 2 ) 2 — Am 2 m 2 \ 

v non = \M = a= K _ P; V IV JL_ ! *1 IZl « (137.02)- 1 , (12) 

Me 2M H m e m p 1 ' V 1 

where the subscripts H, e and p mean the hydrogen, electron and proton, respectively, and fi e is the reduced mass 
of the electron. Solving this equation for the mass of the hydrogen, and it is expanded with respect to the coupling 
constant and m e /m p as 

M H =m p + m e - \o?\ie - \a i fi e + f aVe — - -}-a 6 p e + , (13) 

2 8 8 m p 16 

where the fine structure term (or), which consists of the spin orbit coupling and relativistic corrections, can be derived 
correctly for the ground state. The hyperfine structure term appears less reliable, because the spin-spin interaction is 
not applied to the equation. The equation is also applied to the positronium as 

1 2 1 

—a m P . 

4 64 

Calculating Eq. (|l2) for the deuteron, the result is 0.097 which is the order of the pseudovector coupling constant 
/ 2 /(47r) for the usual definition f/m n = cj-nN n / (2 Ad) JhJ. The ratio of this result to the usual value 0.079 is similar 
to the axial vector coupling constant qa- 

From the above reasonable results, the strong coupling constant can be estimated to a certain extent as 



M posi = 2m e - -a 2 m e - —ot^m e + (14) 



M -- 4 "f = iW (15) 



m q m q 

where ttq has a quark and an anti-quark with the same mass. Since the relativistic velocity is less than 1, the lightest 



quark mass is allowed to be less than 95.4 MeV just as v r q el = \p q \/E q = y\M 2 — 4m 2 |/M 7r o < 1. Hence the coupling 
constant should be estimated from the current quark masses. Now the momentum of the gluon is estimated as 

q 2 = -4p 2 sin 2 (9/2) = -m 2 a 2 sin 2 (9/2), (16) 

in the center of mass frame. Though this equation is defined for scattering processes, if it is applied to QED and 
compared with the Bohr radius, then the relation sin 2 (9/2) = i can be analogized for the ground state (|q| ~< 
V'ioo|l/ r #ioo >= 1/a-Bohr = m e a). 

For a numerical example, the coupling constant a s is estimated as 13.35 for the typical current quark mass m q = 10 
MeV and then the gluon momentum q = 66. 7 'A MeV from Eqs. ([ll]) and (|l6|). The gluon momentum transfered by an 
excited quark can be estimated simply as 201.72 MeV (66.74+134.98) and 836.74 MeV, respectively, at the threshold 
energies of ttq and pq. The running coupling constant a s (q 2 ) — gr^npTA?) can ^ e ex t ra P°l a ted by using unknown 
parameters (3, A for a s (66.74 MeV) = 13.35 and a s (m T ) — 0.35 which is the measured coupling constant at the lowest 
mass scale Ju^pl- The coupling constants can be estimated as a s (201.72 MeV) = 0.988 and a s (836.74 MeV) = 0.451, 
and the expansion parameters of m 2 /q 2 are 0.0025 and 0.00014, respectively, in the energy region of the qq creation. 
Therefore perturbation with respect to the a s and the expansion of m 2 /q 2 are reliable simultaneously in the energy 
region relevant to quark confinement. So, quarks are confined in the hadron by the strong coupling in the non- 
perturbative region and by the linear confinement potential in the perturbative region. If the assumptions of this 
letter hold for the lightest meson, they also hold for all the other hadrons, because the mass of the internal quark 
loop starts from the lightest quark. In QED, on the contrary, the inner most electron of the most heavy atom can 
not transfer a momentum lager than 2m e , because the sum of the binding energy (13.6 x 120 2 eV = 0.196 MeV) for 
the atomic number Z = 120 and the momentum of the photon (0.448 MeV) is estimated as 0.644 MeV. Therefore 
one can not see any confinement potential in QED, because the expansion of m 2 /q 2 is not allowed. 

As a result of the view of this letter, quark masses are not direct observables but parameters to be obtained from 
the hadronic spectra by solving the Dirac equation. So quark masses are unknown exactly. The current quark masses 
support the assumptions of this letter and the mass of the light mesons quite well rather than the constituent quark 
masses. If the quarks are bounded by the usual Coulomb potential, then the mass of the light meson is less than the 
sum of the masses of its constituent two quarks because of their negative kinetic energies. However the quark are 
bounded by the confinement potential of this letter, the mass of the meson is larger than the sum of the masses of 
the two quarks because of their positive kinetic energies. 
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FIG. 1. The Feynman diagrams at tree and one loop levels to contribute to the effective potential. The diagrams for the 
counter terms are not shown. 
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